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We calculate exactly both the microcanonical and canonical thermodynamic functions (TDFs) 
for a one-dimensional model system with piecewise constant Lennard- Jones type pair interactions. 
In the case of an isolated Ai'-particle system, the microcanonical TDFs exhibit (A*" — 1) singular 
(non-analytic) microscopic phase transitions of the formal order N/2, separating A'' energetically 
different evaporation (dissociation) states. In a suitably designed evaporation experiment, these 
types of phase transitions should manifest themselves in the form of pressure and temperature 
oscillations, indicating cooling by evaporation. In the presence of a heat bath (thermostat), such 
oscillations are absent, but the canonical heat capacity shows a characteristic peak, indicating the 
temperature-induced dissociation of the one-dimensional chain. The distribution of complex zeros 
(DOZ) of the canonical partition may be used to identify different degrees of dissociation in the 
canonical ensemble. 

PACS numbers: 05.70.Ce, 05.70.Fh, 36.40.Ei, 64.60. Cn 



INTRODUCTION 

Many macroscopic systems exhibit sudden variations 
of their physical properties (elasticity, conductivity, etc.), 
when one or more external control parameters (e.g^ en- 
ergy E, temperatm'c T) pass certain critical values 
Such phenomena are usually referred to as phase transi- 
tions (PTs). Seminal contributions to the theory of PTs 
are due to, e.g., Mayer et al. 0, 0, IB, Onsager 0, van 



Hove lT I, y a-ng and Lee 
al. llilliilia, Burkhardt 
and Cuesta and Sanchez [18 




Fisher [lOj , Grossmann et 
14, Pettini et al. [ll[ll[l3 

These authors have studied 



in detail canonical ensembles (CEs) and grandcanonical 
ensembles, corresponding to systems in contact with heat 
bath and particle reservoirs. 

Singular, or equivalently non-analytic, PTs are indi- 
cated by a discontinuity in the thermodynamic func- 
tions (TDFs) or one of their derivatives Some- 
times it is also useful to consider smooth PTs, char- 
acterized by a strong but analytic variation in the 
TDFs HO, IH m m. In the presence of a heat bath, 
singular canonical PTs can occur in the thermodynamic 
limit only 0, 0, 0, 0, 0| , whereas finite canonical 
systems may, at best, exhibit smooth PTs |2olEll Ell23t . 
However, the situation changes if the system under con- 
sideration is thermally isolated, corresponding to a mi- 
crocanonical ensemble (MCE). Due to the different phys- 
ical conditions underlying MCE and CE, respectively, 
one can obtain significantly different predictions for sev- 
eral observable quantities [H HI Hi HUll ill. For 
example, in certain cases, microcanonical heat capaci- 
ties can also be negative (e.g. in self-gravitating sys- 
tems) whereas canonical heat capacities are generally 



positive. In particular, as will also be shown below, the 
microcanonical TDFs of finite isolated systems may ex- 
hibit non-analyticities. These singularities reflect evap- 
oration/dissociation phenomena and ma y b e interpreted 
as microscopic PTs in the small system |30j. 

The main objective of this article is to exemplify the 
differences between the MCE and CE and to elucidate 
particular observational consequences by means of a Id- 
model for evaporation. Remarkably, this model sys- 
tem analyzed below allows for calculating exactly both 
the canonical and microcanonical TDFs for an arbitrary 
number of particles. The knowledge of the exact TDFs 
for both ensembles provides the basis for a detailed com- 
parison of observables. Our main results can be summa- 
rized as follows: 

In the case of the MCE the model system exhibits 
{N — 1) singular microscopic PTs, reflected by non- 
analytic kinks in the caloric curve T{E) and the pres- 
sure curve P{E) at certain critical energy values Ek- 
The values Ek can be identified with the binding en- 
ergy of different dissociation states; i.e., the singularities 
(non-analyticities) separate energetically different evap- 
oration phases. These microscopic PTs are accompanied 
by strong temperature oscillations; i.e., the temperature 
of the system decreases when increasing the energy in 
the vicinity of the critical values Ek- This effect cor- 
responds to cooling by evaporation (or dissociation). By 
contrast, non-analytic transitions are absent in the corre- 
sponding CE; i.e., if the system is embedded into a heat 
bath of temperature T. Nevertheless, a smooth PT is 
observed that also persists in the thermodynamic limit - 
even though the existence of a singular macroscopic PT is 
excluded by the (generalized) van Hove theorem 0, E| ■ 
Finally, our study of the distribution of complex zeros 



2 



(DOZ) of the canonical partition function [iolEllElE^ 
suggests that the DOZ may be used to identify different 
degrees of dissociation in the CE. 



THE MODEL 

We consider a one-dimensional model system corre- 
sponding to N identical point particles confined by a 
one-dimensional box of size L. The Hamiltonian reads 



H{p,ci;L,N) = ^ + Ui<i;L,N)^E, (1) 

with q = (gi, . . . , qn) denoting the coordinates and 
p = (pi, . . . ,pn) the conjugate momenta. In the case of 
an isolated system the total energy E is conserved. The 
potential energy U = U-mt + Uhox is determined by the 
interaction potential 

1 ^ 

Uint{q;N) = - U^.Uh-qj\) (2a) 
and the box potential 



f/.„(q;L.A')4°' '^l"'^'" (2b) 

+CXD, otherwise. 



The pair potential is given by 



,0, 



r < dhc 
dhc < r 
r > dhc 



oo, 

pair(?') = { -Uo, dhc <r < dhc + ro, 



(2c) 



where dhc > is the hard-core diameter of a particle with 
respect to pair interactions. The interaction potential 
(|2c|l can be viewed as a simplified Lennard- Jones poten- 
tial. The depth of the potential well is determined by 
the binding energy parameter Uq > and the interaction 
range by the parameter rg, where we shall additionally 
impose that 

< To < dhc- 

The latter condition ensures that particles may inter- 
act with their nearest neighbors only. Furthermore, we 
assume that L > Lmin = {N — l)(dhc + ro), i.e., the vol- 
ume is assumed to be sufhciently large for realizing the 
completely dissociated state, corresponding to J7 = 0. 
The energy E of the system can take values between the 
ground state energy 



Eo = -iN-l)Uo 



and infinity. 



MICROCANONICAL ENSEMBLE 

The microcanonical ensemble (MCE) refers to an iso- 
lated system. Thence, the control parameters are en- 
ergy E, volume L and particle number A'^. The ther- 
modynamic (Hertz) entropy of the MCE is given by 



(3a) 



S{E,L,N) = kBlnn{E,L,N), 
where fc^ is the Boltzmann constant, and 

n{E,L,N)^-^ [ dq/ dvQ{E-H) (3b) 

N ! ft J^N J^N 

the phase volume {h is Planck's constant, and <d{x) = 
for a; < and Q{x) = 1 for a; > 0). Using iV-dimensional 
spherical momentum coordinates, one can rewrite Eq. 
(pbjl as 



C{N) f dq(£;-C/)^/2 e(S-C/), (4a) 

(4b) 



T{N/2)N\Nh^ 

where F denotes the Euler gamma function. For Hamil- 
tonian (P) one can calculate integral (|la|l exactly, yielding 
(see Appendix) 



Af-l 



n^Cj2^k{E + kUof/^ e{E + kUo), 

k=0 

where, for L > {N ~ l)(ro -I- dhc), 



(5a) 



LOkiN, L) 



1=0 



(-1) X 



X [L-{N- l)dhc - 7-0 (TV - 1 - fc + i)] 



N 



(5b) 



Given Eqs. Q, the microcanonical temperature T and 
pressure P are obtained from the standard definitions 
T-i = dS/dE and P/T = dS/dL HHlllllli. For ex- 
ample, for the temperature one finds 



N-l 



^ iOk{E + kUo)'"^Q{E + kUo) 

_f__fe=0 

Y: ujk {E + fcC/o)^/2-i Q{E + fcC/o) 

fc=0 



(6) 



which reduces to the ideal gas law E = NkBT/2 in the 
limit E ^ NUq. It is worthwhile to recall that, for a 
Hamiltonian of the form (Q, the thermal energy © de- 
rived from the Hertz entropy is directly related to the 
microcanonical mean kinetic energy per de gree of free- 
dom by virtue of the equipartition theorem |33l l34| : 



knT 



El 
2m 



= , 2=1,. ..,iV, 



(7) 
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where ( ■ ) denotes the average with respect to the micro- 
canonical probability density function 

Hence, for isolated ergodic systems with an arbitrary par- 
ticle number N, the caloric law T{E) can be measured 
experimentally by monitoring the kinetic energy over a 
sufficiently long time interval (at fixed energy values E) . 

As shown in Fig. ^ the microcanonical caloric law © 
as well as the pressure P{E) exhibit N non-analytic 
points at the energies Ek = —kUo, k = 0, . . . , N — 1, 
separating N energetically different dissociation states 
(all bindin gs i ntact, one binding broken, etc.). The for- 
mal order [SOl Isil Is^ of these non-analyticities equals 
7V/2, i.e., the entropy has continuous derivatives up to 
order {N/2 — 1), but the (A^/2)th derivative becomes 
discontinuous (a similar result was obtained recently by 
Kastner and Schnetz for the mean-field spherical spin 
model 113; see also Gross jl^ for a general discussion 
of differentiability properties of the microcanonical par- 
tition function). Consequently, the 'microscopic (disso- 
ciation) phases' as well as the singularities appear to be 
smoothened out in the thermodynamic limit N —>■ oo. 
Nevertheless, for finite systems - and in particular at 
small densities - the non-analytic behavior is accompa- 
nied by strong variations/oscillations of observable quan- 
tities as temperature and pressure, when continuously 
varying E. Both qualitatively and quantitatively, this 
behavior is analogous to what is usually denoted as a 
'phase transition'. However, since these microscopic non- 
analyticities do not survive in the thermodynamic limit 
(at least for our Id-model), they strictly speaking are not 
covered by the conventional definition of singular macro- 
scopic PTs. We shall, therefore, speak of sing ular (or 
non-analytic) microscopic PTs in the MCE |48l |. 

Let us also briefly discuss the parameter dependence 
of the microcanonical TDFs shown in Fig. As evident 
from Eq. H5a|l . the positions Ek of the singular macro- 
scopic PTs are just proportional to Uq. The amplitude 
of the associated oscillations in T and P does also depend 
on the particle number N and box size L: The strength 
of the oscillations increases for larger values Uq and L, 
but becomes smaller for larger particle numbers N. The 
number and formal order of the PTs, however, only de- 
pend on the particle number N and are independent of 
Uo and L (as long as Uq > and L > Lmin)- Thus, 
qualitatively, the results are independent of the particu- 
lar choice of the model parameters Uq and L. Moreover, 
analogous features can be found in the microcanonical 
caloric curves of Id Lennard-Jones chains [sof . 

It should be mentioned that the exact phase vol- 
ume (|5a|l of our model system resembles in structure 
the phase volume obtained by the Harmonic Superpo- 
sition Method (HSM) applied to Lennard-Jones clusters 



(see, e.g., Doye 01 or Wales and Doye and refer- 
ences therein). The HSM approximates the phase vol- 
ume ^{E) by a sum of ellipsoidal regions around all local 
minima of the potential U lower than the total energy E. 
This method has been successfully applied to describe 
melting phenomena, as e.g. the low-temperature prop- 
erties of 3d Lennard-Jones clusters and their transition 
from a solid-like state, where the cluster only vibrates 
around the ground state configuration, to a liquid-like 
state, where also other locally stable configurations ener- 
getically accessible. The standard HSM approximation, 
however, does not properly account for the contribution 
to the phase volume stemming from (partly) dissociated 
(or gas) states of the cluster, and therefore, is not suit- 
able for describing evaporation phenomena. In particu- 
lar, the HSM does not yield any singular microscopic PTs 
for Id Lennard-Jones chains (where only one locally sta- 
ble configuration, i.e. the ground state, exists) therewith 
contradicting exact analytical and numerical results |3nl | . 
By contrast, the model system discussed here - if consid- 
ered as an approximation to Lennard-Jones chains - does 
reproduce these microscopic PTs related to evaporation 
(but, of course, our model cannot be applied to melting 
processes because it is one-dimensional). 

It is worthwhile to discuss the microscopic PTs and the 
origin of the associated temperature oscillations, as ob- 
served in our model, from a more general point of view. 
Mathematically, microscopic PTs of the above type arise 
whenever the phase volume grows non-smoothly in the 
vicinity of some critical energy value E = Ek- This can 
best be illustrated by considering the energetically ad- 
missible subset of the configuration space 

A{E) ^ {q(^m.^ \Q{E ~U{q)) ^l] . (9) 

The set A{E) consists of all positions space points q = 
(qi, . . . ,qN) that can be occupied by the system at the 
given energy value E. Clearly, the boundary of A, de- 
noted by OA, determines the effective range of the inte- 
gral in Eq. l|Ia|l . Hence, whenever A or dA, respectively, 
change their shape in an irregular (non-analytic) man- 
ner, a non-analyticity in the phase volume O may arise 
(and, hence, in the TDFs). For example, such a irregular 
change in the shape of A occurs when the energy for the 
next dissociation step is crossed, since then some parts 
of the boundary dA suddenly become determined by the 
box potential. 

It remains to be discussed how the temperature oscilla- 
tions - i.e., the regions with negative heat capacit y (a lso 
known as 'S-bends' or van der Waals-type loops |42j) - 
arise: In the vicinity of the dissociation energy Ek, the set 
A and, thus, also O and S grow very rapidly, thereby giv- 
ing rise to a drop-off in temperature. Geometrically, this 
can be viewed as a sudden increase of the 'effective di- 
mensionality' of yi. Here, 'effective dimensionality' refers 
to the number of orthogonal configuration space direc- 
tions in which A has an extent comparable to the system 
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size L. Hence, typically, the temperature oscillations ap- 
pear more pronounced for larger values of L. From the 
physical point of view, the temperature decrease after the 
fcth dissociation step just means that for energy values 
slightly larger than Ek the dissociated fragments have 
very little kinetic energy (since most of the energy has 
already been used to break the binding). The larger 
the system the less likely it is that the fragments tem- 
porarily recombine in a state of high kinetic energy; i.e., 
from a probabilistic standpoint, the average lO is then 
dominated by phase space regions of low kinetic energy. 
With regard to practical applications, this means that 
one could cool such a small isolated system of bound par- 
ticles by injecting energy until the fragmentation process 
sets in (cooling by evaporation/dissociation). 

The above described features of microscopic PTs are 
generic and shared by all physical systems that exhibit 
dissociation and evaporation (e.g. similar microscopic 
PTs and temperature oscillations are also found for small 
Id-Lennard- Jones molecules [sof'). In particular, micro- 
scopic PTs should become more pronounced in two or 
three dimensions (since then effective dimension of A 
grows even more rapidly at the dissociation levels) and 
also be observable in quantum systems. By virtue of dis- 
sociation experiments with small particle numbers and 
very low densities (similar to those of Schmidt et al. j 
but without heat bath), one should, in principle, be able 
to detect the oscillating behavior e.g. in temperature and 
pressure curves. However, to actually observe such oscil- 
lations one has to realize the requirements of the MCE, 
i.e. a thermally isolated system with regulated energy 
injection. Furthermore, due to the microscopic origin of 
the oscillations and the requirement of a relatively low 
particle density, a high sensitivity of the velocity and 
force measurements and long measuring time spans will 
be necessary. 
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FIG. 1: (a) Microcanonical temperature T and (b) pres- 
sure P as a function of energy per particle v — E/N for 
a (reduced) density n = N/[L - {N - l)dhc] = 0.001/ro and 
different number of particles N — 5 (dashed line), TV = 15 
(dotted), and A'^ = 500 (solid). Note that each of the curves 
is (A'^/2 — 2) times differentiable. 



with LOk (-/V, L) given by Eq. I)5b|l . Mean energy and pres- 
sure of the CE are defined by E = —d{lnZc)/d(3 and 
P = —dF/dL, yielding e.g. the canonical caloric law 



CANONICAL ENSEMBLE 

Employing the canonical ensemble (CE) is appropriate, 
if the system under consideration is in thermal contact 
with a much larger system (heat bath), as e.g. realized 
in dusty cluster experiments |42l |. The relevant thermo- 
dynamic potential is the free energy 



F(/3,L,iV) = -/3-ilnZc(/5,L,7V), 



(10) 



where Zc is the canonical partition function. The ex- 
ternal control variables are now the inverse temperature 
(3 = {kBT)~^ of the heat bath, the volume L and the par- 
ticle number N . For the above model, Zc can be exactly 
calculated, analogous to Eq. ©, as 



1 / 27rm 



N\ \ (3h 



N/2 N-l 

Y^oj.eP^^- 

fc=0 



(11) 



E 



N 



E^Ao^keP'^Uo 



(12) 



Figures H (a) and (b) show E{T) and P{T) 
for different values of the reduced particle density 
n = N/[L — {N — l)c?hc]- In contrast to the microcanon- 
ical pressure [Fig. ^ (b)], the canonical pressure is a 
monotonous function of T or E, respectively. In the ther- 
modynamic limit, microcanonical and canonical caloric 
curves become nearly indistinguishable. The canoni- 
cal heat capacity cl = BE/dT exhibits a strong (non- 
singular) peak in the temperature region, where dissoci- 
ation occurs [Fig. 121(c)]. If observed in an experimentally 
measured curve, such behavior would possibly be inter- 
preted as a PT. For decreasing density n, the position 
of the maximum of cl moves closer to T = 0, while its 
height increases rapidly. Furthermore, our results indi- 
cate that for iV > 15 the TDFs v and become virtually 
independent of N . The (non-singular) peak in the heat 
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FIG. 2: (a) Canonical mean energy per particle u — E/N , (b) 
pressure P and (c) specific heat capacity cl ~ dE/dT (loga- 
rithmic scale) as a function of temperature T for A'^ = 15 par- 
ticles and different values of the reduced density n — 
(dashed line), n = 10~^/ro (dotted), and n = 10"*^ /ro (full). 



capacity persists in the thermodynamic limit (analogous 
to the Id Ising model [i^). 

To obtain a more detailed characterization of the disso- 
ciation process in the CE, we next study the distribution 
of complex zeros (DOZ) of Zq- As evident from Eq. 
the only relevant configurational part of Zq (/3) is a poly- 
nomial of degree {N — 1) in z = e^^°, and, therefore, has 
{N — 1) complex Fisher zeros P^r branch of the loga- 
rithm. This quasi-polynomial structure is a consequence 
of the fact that, for our specific model, the configuration 



space [0, L]^ can be partitioned into regions of equal to- 
tal binding energy Ek = —kUo, k = 0, . . . , N — 1 (see 
Appendix) . Since all zeros can be obtained by adding in- 
teger multiples of 27ri to the zeros (3k of the main branch 
(for which 9 (Pk) = tt), it suffices to discuss to the main 
branch only while bearing in mind that each zero (3k in 
the main branch is associated with an infinite set of zeros 
{(3k + 2Tii s\s S Z}. Ordering the zeros according to their 
real parts, 5ft(/3o) < • . . < 'SI{(3n-i), we find that the re- 
gion of the Ci-peak is well-described by the temperature 
interval [3?(/3Ar_i)-i, 5R(/3o)~^]. 

The asymptotic behavior of the DOZ for N ^ oo 
may be used to characterize the parameters (critical 
ternperature, order, etc.) of singular macroscopic PTs 
[ilaElIilElEa- in our model, the Fisher zeros are 
located at least a distance it away from the the real (3- 
axis regardless of particle number N (see Fig.OJ. Hence, 
the zeros cannot converge to a point on the real /3-axis. 
The peculiar position of the zeros thus ensures agreement 
with the (generalized) van Hove theorem 0,^3 which ex- 
cludes the existence of a singular macroscopic PT in our 
model. 

The DOZ has also been employed to study their finite 
size analogs of macroscopic PTs m m . Apply- 
ing these methods to our model, one may interpret the 
smooth phase transition observed in the CE as a super- 
position of (TV — 1) smooth microscopic 'first-order' phase 
transitions indicated by the (A^ — 1) sets of zeros {(3k + 
2-Kis\s £ Z}. One can then use {lR{(3o), . . . ,^{(3n-i)} 
to distinguish (define) different dissociation states in the 
CE. Thus, our results suggest that the DOZ encodes 
detailed information about the observed smooth phase 
transition and the energetically different degrees of dis- 
sociation even if there is no singular PT in the thermo- 
dynamic limit. In particular, since microcanonical phase 
volume and canonical partition function can be mapped 
onto each other via the Laplace transformation |4J, |45| , 
one may speculate that there exists a direct mathematical 
link between the singular microscopic PTs in the MCE 
and the DOZ in the CE. 



SUMMARY 

We have studied a simple Id-evaporation model with 
nearest neighbor interaction potentials characterized by 
a hard-core repulsive part and piecewise constant short- 
range attraction. By analyzing the exact TDFs of this 
model, it was shown that, in the case of a thermally iso- 
lated system (MCE), the microcanonical caloric and pres- 
sure laws exhibit singularities, separating different disso- 
ciation states. The formal order of these non-analyticities 
increases as the particle number increases. Hence, the 
microscopic PTs vanish in the thermodynamic limit and 
are intrinsically different from the topqlogically induced 
macroscopic PTs discussed in Refs. M hj, hj, UO, k7| . 
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FIG. 3: DOZ of the canonical partition function Zc(/3) for 
N — 5 (open symbols) and = 15 (filled symbols) and 
n — 10~^/ro (triangles), and n = (diamonds). For 

better visibility, for N — 5 (15) only the first few branches 
of zeros with positive (negative) imaginary part are shown 
(further branches can be obtained by shifting with multiples 
of 2ni in the vertical direction; generally, the Fisher zeros are 
symmetric with respect to the real axis 12^). 



For sufficiently low particle numbers and densities, the 
microscopic PTs are accompanied by strong oscillations 
of the temperature (mean kinetic energy) and pressure. 
These oscillations arise from a rapid change of the phase 
volume near the dissociation thresholds. They are a 
generic feature of particle systems with Lennard-Jones 
like interaction potentials [sol l39| . They are not re- 
stricted to one space dimension, but expected to be even 
stronger in two and three space dimensions, and should 
also be observable in quantum systems. In a suitably 
designed dissociation experiment, one should therefore, 
in principle, be able to detect the oscillating behavior 
e.g. in temperature and pressure curves. In particular, 
such temperature oscillations may provide the possibil- 
ity to cool a small isolated system by means of regulated 
energy injection (cooling by evaporation). 

If the model system is coupled to a heat bath (CE), a 
smooth PT is observed, but no singularities are found 
0, Nevertheless, the DOZ seemingly permits to 

quantify the temperature range of the PT and to define 
different dissociation states. 

We thus conclude this paper with two important ques- 
tions, which need to be answered in the future: Is it 
possible to design an experiment that allows to observe 
singular microscopic phase transitions in finite size sys- 
tems (or, at least, oscillatory behavior of thermodynamic 
observables such as pressure) , as predicted by the micro- 
canonical statistical theory? Can one find a direct mathe- 
matical link between singularities in microcanonical par- 
tition function of finite systems and the DOZ of the cor- 
responding canonical partition function - and, thus, be- 
tween microscopic and macroscopic phase transitions in 
arbitrary space dimensions? 
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For k 



0, . . . , N — 1 the family of sets 



Sfe = {x e [0, A]^ I C/'(x) = -kUo} constitutes a 
partition (disjoint cover) of the configuration space 
permitted by the transformed box volume [0, A]. Thus 
we can rewrite Eq. (|4a|) in the form 

n = c{N) + kUof^ e{E + kUo) voi(gfc). 

In order to calculate vol(S/c) we first note that 

Uk=yoli9k) = iV!.vol(g+), 
where (writing |xg[0,A]^| ...} = {...} from now on) 

St = {xi < X2 < ■■ . < XN A C/(q) = -kUo} . 
Next we note that 



where 

voi(g++ 



{xi < X2 < ■ ■ ■ < xn} n 

{{x2 - Xi < To) A . . . A {xk+1 ~Xk < tq) a 
{xk+2 - Xk+1 >ro) A... A (xN - xjv-i > ro)}, 

dxN / da;Ar_i x • ■ • 

{N-l-k)ra J (N ~2-k)ro 

Xk+i-ra f-Xk + 2-ra 

dxk+2 / dxk+1 X 

ro "'0 
Xk+1 r^a /•X2 

dxk ■ ■ ■ dx2 / dxi. 

Xk+i-ro]+ Jlx3-ro]+ Jlx2-ro] + 



The positive part [x]+ = max{0, a;} satisfies 

.]n+l 



dx \x]l 



X 



n + 1' 

la; — c| , — c|_i_ = \x ~ 2c] + , 



n eN; 
c> 0. 



(14a) 
(14b) 



J+ - H + 

By virtue of these identities, we can write Eq. (|I4|I as 



APPENDIX 



To calculate integral (|4a|l . we first eliminate the hard- 
core part of the interaction potential by virtue of the 
transformation L ^ X — L — (N — l)dhc, li '—^ Xi — 
qi — riidhc G [0, A], here rii equals the number of particles j 
with Qj < qi. With these definitions the potential energy 
can be rewritten as: 

C/'(x; \,N) = \Y^U'^^,{\x, - x,\) + C/box(x; A,iV) , 
where the sum ^ goes over nearest neighbors only, and 



-Uo, 0<r<ro, 
0, r > ro . 



vol(Sfc ^) = / dxN / da;jv-i x • • • 

J(N~l-k}ro J{N-2-k)ro 
fXk+s-ro rXk + 2-rQ 

dxk+2 I dxk+1 Kk, 

iro Jo 



where 



i=0 



Performing the remaining integrations and reversing the 
transformation _L i— > A, one obtains the final result Ij5b|l . 



